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We present a class of thermodynamic systems with constant thermodynamic curvature which, 
within the context of geometric approaches of thermodynamics, can be interpreted as constant 
thermodynamic interaction among their components. In particular, for systems constrained by the 
vanishing of the Hessian curvature we write down the systems of partial differential equations. In 
such a case it is possible to find a subset of solutions lying on a circumference in an abstract space 
constructed from the first derivatives of the isothermal coordinates. We conjecture that solutions 
on the characteristic circumference are of physical relevance, separating them from those of pure 
mathematical interest. We present the case of a one-parameter family of fundamental relations that 
- when lying in the circumference - describe a polytropic fluid. 
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I. INTRODUCTION 



The study of physical systems that admit a geometric description in terms of Riemannian manifolds is an interesting 
and timely subject. Over the last few years, there has been a number of efforts towards an integrated description 
of thermodynamics in terms of Legendre invariant quantities. In particular, analogously to the case of field theories, 
it has been argued that the curvature of the appropriate manifold should be linked to the notion of thermodynamic 
interaction [ij. There has been numerous proposals in this direction. On the one hand, there are the conformally 
related metric theories of Ruppeiner and Weinhold , where the metric takes the form of a Hessian of the extensive 
parameters in the entropy and energy representations, respectively 3i, 3:]. However, both fail to comply with the spirit 
of the geometric construction of field theories, i.e. those are not invariant under the natural set of transformations in 
thermodynamics. On the other hand, the Geometrothermodynamics programme (GTD) has successfully managed to 
provide us with a set of metrics which are independent of the potential used [ij and the fundamental representation 
one uses [1, @ . 

In the GTD programme one posits that the physical information about a thermodynamic system cannot depend 
on the potential used to describe it and that such information is encoded in the curvature of the maximal integral 
manifold of the Pfaffian system defining the first law of thermodynamics [c.f. equation ([T])]. We call such a manifold 
the space of equilibrium states. The curvature of such manifold is obtained from the first fundamental form, induced 
but a Legendre invariant metric of the contact manifold specified by the Pfaffian form. In section|lT]we present a brief 
review of the programme, in particular, we centre our attention on a metric whose curvature does not depend on the 
fundamental representation. 

Thus far, the GTD formalism has been applied to a number of thermodynamic systems in order to test the 
consistency of the programme. In this work, we take a step forward and set ourselves to the task of finding those 
systems which exhibit constant thermodynamic interaction. That is, we will find the class of fundamental functions 
producing a manifold of constant curvature. We will restrict ourselves to the two-dimensional case, i.e. we consider 
only systems with two degrees of freedom. These systems are interesting from the mathematical point of view since 
the any two dimensional metric can be cast into a conformally fiat form. 

The paper is organised as follows. In section II we present a brief account of the Geometrothermodynamics 
programme and two-dimensional thermodynamic systems of constant curvature. In section HI we analyse the system 
of partial differential equations to find the set of isothermal coordinates for metrics with vanishing Hessian curvature. 
There we propose a criterion to single out physical fundamental relations based on a circumference-like equation in an 
abstract space related to the system of differential equations for the isothermal coordinates. To close this section we 
present an example illustrating these matters. Finally, in section IV we present a summary of the results and address 
future work on the subject. 



II. TWO-DIMENSIONAL THERMODYNAMIC SYSTEMS OF CONSTANT CURVATURE 



The GTD programme promotes the natural formalism of thermodynamics in terms of contact manifolds to a 
Legendre invariant Riemannian structure. Let us begin with a brief review of the programme by considering the 
case of two thermodynamic degrees of freedom. In this case, we need a five dimensional manifold which admits a 
set of local coordinates corresponding to the collection of extensive and intensive variables - denoted by and pi, 
respectively - together with the thermodynamic potential, such that the kernel of the 1-form 

e = d$ -pidqi -p2dg2, (1) 

generates a maximally non-integrable set of hyperplanes, f C TT- A manifold T together with the 1-form 9 is called 
a contact manifold. In the present case we refer to it as the thermodynamic phase-space. Of special interest is the 
maximal integral sub-manifold, £ C T, i.e. the largest sub-manifold which can be embedded in 7" such that T£ C ^. 
It is easy to see that this is a two-dimensional manifold defined by the first law of thermodynamics 

= pidqi + p2dq2, where $ = $(91,92) and pi = —— = ^,. (2) 

oqi 

Thus, we see that if we know the fundamental function $ — $((71, 92), then we know how £ is embedded in T- We 
call the sub-manifold £ the space of equilibrium states. 

In addition, the GTD programme introduces a metric structure for the thermodynamic phase space. Such a structure 
is constructed in order to satisfy the criterion of Legendre invariance, i.e. Legendre transformations correspond to 
isometrics. Within the GTD programme there have been two distinct classes of metrics which have been studied 
according to their invariance properties, those which are invariant under every possible Legendre transformation and 
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those which are only invariant under total Legendre transformations. The metric structure of T induces a Riemannian 
metric on £, its first fundamental form, whose intrinsic curvature is associated with the thermodynamic interaction 
of the system. In our two-dimensional scenario, this whole information is contained in the curvature scalar of £. 

If the curvature of the space of equilibrium states is to give a faithful account of the thermodynamic interaction, 
it should not depend on the choice of fundamental representation, i.e. one is free to work in the energy or entropy 
representation indistinctly. It has been shown that the metric compatible with both, Legendre and representation 
invariance is 

G'' =9(8)6 + ^— (dgi® dpi +dg2«)dp2). (3) 

92^2 

Thus, the induced metric on £ is simply given by 

g^=niq,,q2)h. (4) 

Here h is the Hessian metric 

h = dqi (g) dqi + ($42 + <l',2i) dqi ® dq2 + $,22 <iq2 (g) dq2, (5) 

where we have used a coma to denote partial differentiation with respect to the corresponding coordinate function of 
£ and the conformal factor is given by 

n{qi,q2) = ^. (6) 

The interested reader in the derivation of the metric (j4]) is referred to 4] and to 5] for applications to ordinary 
thermodynamic systems. 

Note that the components of the metric (jlj depend on the second derivatives of the fundamental function $ but 
are otherwise unspecified. It is an interesting exercise to find a class of fundamental functions for which the space of 
equilibrium states £ becomes an Einstein manifold for the metric (|4]). That is, we look for solutions of the system 

R\t-Kg\„ (7) 

where R\ij is the Ricci tensor associated with and X is a constant, which in the present case corresponds to the 
Gaussian curvature of £. 

It is worth noting that equation ([7]) represents a system of three, third order, non-linear partial differential equations 
for the thermodynamic potential $. Indeed, it is straightforward to show that in two dimensions, the fourth order 
terms in the curvature exactly cancel whenever the metric is the Hessian of a scalar function. 

We can reduce the system ^ by raising one of the indexes to obtain 

RI'^kV. (8) 

Thus, the system reduces to the single PDE 

-4i^ 7^5/^ r with i,j-l,2. (9) 

f^l-'[qi,q2) 

Here p is the determinant of the metric ^ given by the expression 

p = n\qi,q2) {^,11^,22-^^2), (10) 

andthe Ihs of is 

= <^>% {A $2^1 + B $,11 - 2g2$,2ii*a2 + C $42 + D $^22) 

+ *,2 (g2^,22«',222*41 + E $41 + 2g2$3^2*,221 - g2*42*,22$,21l) 

- 2qi p {^%2'^M - 'i>'l2'i>,22) , (11) 

A= - 52^,222 - 2$^22, (12) 

B ^ 2$2^2 + 392$,12$,221 - g2«',211$,22 + g|$^221 ^ 9|5',222«',211 , (13) 

C= -g2*,221*,211+g2*.lll*,222 +g2«',22«',lll, (14) 

= -92*,111*,221 +g|«'211, (15) 

E ^ - 2g|$,22*,12^,221 - <7i*^12*,222 + <7i^'^22^,211 • (16) 



where 
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Motivated by the results of a previous work by the authors (c.f. Section III.D in Q and @), we know that a 
solution to is given by the fundamental relation 

$ = <f>olog(<z?+cg«). (17) 

Since we are working in two dimensions, the Gaussian and scalar curvature are proportional and we see that the 
constancy of K is satisfied and has the value 

4 a — 1 

Thus, we can propose a general solution of the form 

$ = /(e'?i+X'?2). (19) 

Here / is a sufficiently differentiable function of the sum of the extensive parameters where ^ and x ^^re constants. 
This type of ansatz does solve (dJ). However, a quick inspection to the metric determinant reveals the degeneracy of 
this case [c.f. equation ([2T|) . below]. Therefore, let us propose the more general solution 

^^f{Cqi+Xq2), (20) 

where a is a constant different from one. Now, the metric determinant is in general different from zero and has the 
form 

?+~\,, [/'+e«g?(eg?+xg2)r], (21) 

X 92 iff 

where /' and /" are the first and second total derivatives of the fundamental relation ((20)) evaluated at {£, + X Q2) ■ 
Now we can clearly see the degeneracy for a = 1. 

Substituting our ansatz ([20)) into the equation ([8]) we obtain again the same result as in the case of (flT)) . i.e 
the Gaussian curvature is the same constant, equation ()18p . regardless of the particular form of the function / as 
long as the argument is (^ 9" +X (Zf)- Therefore, the generalised Chaplygin gas, equation (fT7)) belongs to a class 
of thermodynamic systems with the same type of interaction given by (|20p. Moreover, the Hessian metric for the 
logarithmic form of this type of fundamental relation has vanishing curvature. In this case it becomes a simpler 
problem to find the set of isothermal coordinates for the space of equilibrium states. 



III. ISOTHERMAL COORDINATES 



It is a well known result that every two-dimensional Riemannian manifold is conformally flat. That is, we can 
always find a set of coordinates for which the metric takes the form 

g = VL'^{x,y)g' , where g' = dx ® dx + dy ® dy. (22) 

Such a coordinate system is called isothermal. In this section we find the isothermal coordinates for the space of 
equilibrium states {£,g^) under the assumption of the Hessian flatness, i.e. by demanding that the curvature scalar 
of the Hessian part of the metric (|4]) vanishes. 

Let us consider the diffeomorphism if : £ ^ £ accounting for the change of coordinates x — a; (91,(72) and y — 
2/(911 92)- Then we can pull-back the metric [c.f. equation ()22p ] and solve the equation 

ip*g^ -h = (23) 



for the coordinate functions x, y and the thermodynamic potential $. This will provide us with a family of ther- 
modynamic fundamental relations with zero Hessian curvature together with their isothermal coordinates. Equation 
(|^ above corresponds to the system of equations 

(24) 

2^,1^,2 + y,iy,2 = $,12, (25) 
4 + y% = $,22. (26) 
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There is a large family of solutions for such a system. In particular, for separable fundamental relations, 

<^ = S{qi)+T{q2), 
we have that the change of coordinates is given by 



V(l - c2) S" dqi + c 



y 



(1-C2)f d(72-C 



V T dq2, 



where c is a constant and the primes and dots denote differentiation with respect to qi and (72, respectively. 
A more general solution is given by those fundamental functions satisfying the third order system of PDE's 



5>.211 = 



$,222*^12 



and $,1,22 — 



*&,222*l'l,2 



22 



In this case, the isothermal coordinates must satisfy the system 
1 



^'1 



(2$, 122^,12^,2 + ^,11*I',22 — 2; 2*1', 11 



1 $.2222:, 2 



2; 22 = — - 

' 2 $,22 



and 



$.22 - x\dq2 + - 



22 



22:, 22:, 12 - 'I>,221 



dq2 - 2x,ix,2 + 2$42 



dqi 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 

(33) 



One can verify that a fundamental relation of the form ([T9| is a solution of ([30| . We have seen that this type of 
functions generates degenerate Hessian metrics. However, we can use them to learn some properties about the space 
of solutions of the system ((24)) - ([26]) . For example, consider the fundamental relations given by 



$ = log(C(7i +XQ2)- 

In this case we can solve the pair of equations for x, i.e. equations ()31|) and p2p to obtain 



cc = clog ( 92 + -91 
X 



and substitution in ([55]) yields 



y 



v/-(l + c2)log(e(7i +X(72). 



(34) 
(35) 
(36) 



Thus we see that, indeed, this type of fundamental relation fails to produce a real change of coordinates satisfying 
(1^51) . Moreover, note that we can find particular solutions to the system - if we restrict ourselves to a 
circumference in an abstract XY plane. Thus, we have 



where 



2;.2j 



(37) 

(38) 
(39) 
(40) 

From this point of view, we observe that the fundamental relation ([34]) corresponds to an 'imgainary' radius of the 
circumference (l37l). that is 



=(.T, 

Y'={yj+y,2f 

=$41 + 2$42 + *,22- 



{i + xY 
{ill +xq2Y 



(41) 



This is not surprising since we knew that the Hessian corresponding to this fundamental relation is degenerate and 
thus the system is not well posed except for the case ^ = — x = 1, for which i? = 0. 

We can use this geometric construction to probe the space of solutions for a fixed fundamental relation by noting 
that a solution to the system (IM)) - (pS)) must lie on the circumference associated with the particular fundamental 
relation we use [c.f. equation (|41l) ]. but not every solution lying on the circumference solves the system we are probing. 
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A. Example 

To see how this construction works, let us chose a family of fundamental relations in the form of equation p7)) 
parametrised by the exponent a. We work in the entropy representation using molar quantities. Thus we set by 
Qi = u the specific energy and q2 = v is the specific volume. The fundamental relation is written as 

s„ = logK + i;"). (42) 

Each of these functions defines a Hessian metric of zero curvature and a natural metric of constant thermodynamic 
interaction [c.f. equation (1181) ]. The change to isothermal coordinates for this type of functions cannot be expressed 
analytically. However, we can use the circumference to classify the various types of differential equations obtained for 
each value of a. 

The squared radii of the circles associated with each function is given by 

Rl = -—^ ru2a„2 _ ^a^a _ ^^^2 _ ^aUV + {a - 1)) + V^"u^] (43) 

Thus, fixing the volume of the system, we find the subset of fundamental functions for which the PDE represented 
by the circumference p7p is well defined. The three qualitatively different types of behaviour are depicted in figure 
[TJ It is easy to observe that for a > the PDE is ill defined, thus, only fundamental relations for which a < are 
meaningful. This corresponds to a positive thermodynamic curvature as can be seen from (jl8p . Furthermore, noticing 
the symmetry of u and v in the expression for the circumference radius p3)) . it is easy to observe that the PDE also 
restrict the domain of the thermodynamic variables to u,v > 0. The fundamental relation with a < describes a 
polytropic fluid with equation of state given by 

P = pi-", (44) 

where p = ^ is the energy density of the fluid. It is a simple task to obtain the heat capacity at constant volume for 
these systems 

= , . (45) 

From this expression we observe that the heat capacity remains finite for any value of the thermodynamic variables 
and is always negative whenever a < 0. 



IV. CLOSING REMARKS 



In this paper we studied two-dimensional Einstein manifolds for the Geometrothermodynamics programme. We 
found the differential equation that must be satisfied by the fundamental relation in order to describe a system 
with constant thermodynamic interacion, i.e. a fundamental relation producing a representation invariant metric 
whose associated curvature is constant. In particular, building on previous work (c.f. reference Q) we analysed the 
one-parameter family of fundamental relations given by equation ()20p . 

Noting the conformal structure of [see equation Q], we centre our study in the class of functions whose associated 
Hessian metric has vanishing curvature. With this assumption, we set up the system of differential equations defining 
the set of isothermal coordinates. As expected, we found an analytic expression for the change of coordinates for the 
case in which the fundamental function is separable. An interesting exercise allowed us to explore some properties of 
the space of solutions of i.e. the set of isothermal coordinates with their corresponding fundamental relation of 
vanishing Hessian curvature. 

Observing the algebraic structure of the system (IMl) - (PE)) . we note that there will be a class of fundamental 
relations satisfying the Hessian curvature constraint for which we can build a characteristic circumference on which 
the solutions of the PDE system lie. This can be done whenever the derivatives of the fundamental relation define a 
positive squared radial function. Moreover, we conjecture that only the set of fundamental relations for which such a 
construction is possible can describe physical systems of constant thermodynamic interaction. In particular, we work 
out the example given by (jl7p . Here, we work in specific thermodynamic variables in the entropy representation of 
a system characterised by the exponent a. Indeed, only those systems for which a < correspond to a polytropic 
fundamental relation. 

In sum, we have analysed a particular class of fundamental relations of constant thermodynamic curvature. It 
remains to explore the larger class of functions within the set of solutions of equation (1231) and study their thermody- 
namic implications. This will be done in a forthcoming article. 



FIG. 1. The three quahtatively different types of radial functions. The horizontal axis corresponds to the energy, while the 

vertical represents the radius of the circumference. Wc observe that only the curves corresponding to negative values of the 
exponent a are positive definite in the full physical domain, i.e. positive values of energy and volume. 




FIG. 2. These plots correspond to the radial function for — 1 < a < 1. We observe that only those values of a < can be 
associated with physical systems. 



8 



ACKNOWLEDGEMENTS 

ACGP is funded by a TWAS-CONACYT postdoctoral grant. CSLM is thankful to CONACYT, postdoctoral Grant 
No. 290679_UNAM. FN acknowledges support from DGAPA-UNAM (postdoctoral fellowship) 



[1] H. Qucvcdo, Geometrotherrnodynamics, J. Math. Phys. 48, 13506 (2007). 

[2] G. Ruppeiner, Thermodynamics: A Riemannian geometric model, Phys. Rev. A 20, 1608 (1979). 
[3] F. Weinhold, Metric geometry of equilibrium thermodynamics, J. Chem. Phys. 63, 2479 (1975). 

[4] A. Bravetti, C. S. Lopez-Monsalvo, F. Nettel and H. Quevedo, The conformal metric structure of geometrothermodynamics, 
J. Math. Phys. 54, In press (2013); arXiv:1302.6928 [math-ph]. 

[5] H. Qucvcdo, F. Ncttcl, C. S. Lopcz-Monsalvo and A. Bravetti, Representation invariant Geometrotherrnodynamics: appli- 
cations to ordinary thermodynamic systems, (2013); arXiv: 1303. 1428 [math-ph]. 

[6] A. Aviles, A. BasterrechearAlmodovar, L. Campuzano and H. Quevedo, Extending the generalized Chaplygin gas model by 
using geometrotherrnodynamics, Phys. Rev. D 86, 063508 (2012). 



